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We propose a dielectrophoresis model for phase-separated manganites. Without increase of the 
fraction of metallic phase, an insulator-metal transition occurs when a uniform electric field applied 
across the system exceeds a threshold value. Driven by the dielectrophoretic force, the metallic clus- 
ters reconfigure themselves into stripes along the direction of electric field, leading to the filamentous 
percolation. This process, which is time-dependent, irreversible and anisotropic, is a probable origin 
of the colossal electroresistance in manganites. 
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Introduction. Phase separation (PS) can be found in 
various condensed matter systems. In the past half a 
century, the effects of electric field on phase-separated 
liquids have been cross-fertilized by physicists, chemists 
and biologists Jj, i3i; iSj] . Among them, a distinct phe- 
nomenon is dielectrophoresis : the mixing or demixing 
of neutral materials with different dielectric constants un- 
der a non-uniform electric field [2j l3|j l4, iSi] • It has also 
found wide application in microbiology, nanotechnology, 
polymer processing, etc. 0)0]- However, dielectrophore- 
sis would seem unlikely in solids, which normally lack the 
necessary condition of flowability. 

As for solids, PS has been demonstrated to be crucial 
to understanding various exotic phenomena in strongly 
correlated electronic materials [Sj, manganites being a 
typical example. In manganites, PS generally consists 
of ferromagnetic (FM) metallic phase and insulating 
phase which is usually antiferromagnetic and charge- 
ordered (CO) 0]. One can undertake to understand 
the colossal electroresistance (CER), first observed in 
Pro.yCao.sMnOs where an applied electric field 

could reduce the resistance enormously. Subsequent ex- 
periments found a simultaneous uprush of magnetization, 
implying that the CER might be caused by the collapse of 
CO phase to FM phase [llj. Ahhough this CO-FM col- 
lapse scenario can lead to the CER, the estimated thresh- 
old of electric field is one order of magnitude larger than 
the experiment V^. Moreover, in a recent intriguing 



experiment on Lao.225Pro 4Cao.375Mn03, no appreciable 
change of the magnetization was detected when the elec- 
tric field induced a resistance drop of as large as four 
orders of magnitude |13|]. This anomalous CER calls for 
a further understanding of PS beyond the equilibrium 
ground state in manganites. 

In this Letter, we propose a dielectrophoresis model to 
understand the PS dynamics associated with the CER 
in manganites. At the beginning, it is essential to 
demonstrate the possibility of dielectrophoresis in man- 



ganites being solid. Unlike liquid mixtures, the differ- 
ent phases in a prototype phase-separated manganite 
have the same chemical composition. The migration 
of eg electrons rather than cations makes the different 
phases in manganites 'flowable' in a wide temperature 
range [li, 15, la|, satisfying the first condition of dielec- 
trophoresis: ambulatory phases. The second condition, 
non-uniform electric field, can be solely self-generated 
by the PS, as shown in Fig. 1, without external con- 
trol which is necessary for dielectrophoresis in liquids 0] . 
Therefore, dielectrophoresis is not impossible in mangan- 
ites although they are solids. 

The model. As shown in Fig. 1, our model is based on 
a two-dimensional square lattice (L x L, L = 64) with 
a voltage drop V applied in the y direction. A fraction 
Pm of the lattice sites are metallic, the others insulating. 
This binary mixture picture is adequate to describe the 
large scale PS consisting of FM and CO phases A 
standard Metropolis algorithm is employed in our Monte 
Carlo simulation with the dielectrophoresis mechanism: 
exchange between nearest neighboring (NN) sites with 
the probability determined by the free energy change. 
Consideringthe free energy of electric field Fe in the 
medium [J13|> the total free energy F in our model can 
be written as, 

F = Fe + Fs = -^J e{r)E\r)dr + j^AdS, (1) 

where e(r) denotes the dielectric constant and E{r) is 
the local field at point r. e(r) equals Em (e/) if the phase 
at r is metallic (insulating), where Em (e/) is the dielec- 
tric constant of the metallic (insulating) phase. In real 
manganites, many factors can play the role of a barrier 
in the transition between metastable states, e.g. domain 
interfaces and defects 1^, 1^, and here we simplify the 
physics by characterizing them with the interface energy 
Fs [ul in Eq. (1). Fs is proportional to the interface 
area S between metallic and insulating regions with a 



2 




rivj. i: ^i^oior ^^mmej ijeiu: me pnaae separaueu muuice con- f 

sisted by metallic (red) and insulating (cyan) sites. Right: 0.00 0.02 0.04 0.06 0.08 0.10 



the corresponding color image map of V{r). E 



coefficient A. Within the system, E(r) is obtained us- 
ing the resistor- network (RN) model 3, 2^. In the RN 
model, three types of resistors are distributed according 
to the links between NN sites: Rm between metallic sites, 
Ri between insulating sites, and Rmi between a metallic 
and an insulating site [lOl, set as {Rm + Ri)/'^ here. The 
voltage of each site V{r) in the RN can be exactly solved 
by KirchhofF equations, as shown in Fig. 1. 

Our simulation begins with a zero-field quench prepro- 
cess from a high temperature Tq to T^. Then at fixed, 
the voltage V changes between zero and a maximum 
value Vm, following an increasing-decreasing-increasing 
(IDI) path with a fixed rate (characterized by Monte 
Carlo steps (MCS)). Given the strong contrast between 
the electric properties of FM and CO phases 9, 2l|, both 
the ratios Ri / Rm and em/si are quite large. Their val- 
ues strongly depend on material, as well as doping and 
temperature. Experimentally, they can be determined, 
by comparing the metal-rich limit and insulator-rich limit 
@) El. To accelerate our simulation, a special case: 
/ = Ri/Rm = £ai/£i (/ ^ 1) will be used. More general 
cases without such restriction have also been tested. It 
is found that the physical picture is unaffected as long as 
/ ^ 1. The values of parameters used in our simulations 
are listed in Table. I, if not noted explicitly. 

Results. The resistance Ry along the direction of the 
electric field is shown in Fig. 2 as a function of the ex- 
ternal electric field E {E ^ V/{L - 1)). Initially, the 
resistance is only weakly dependent on E. Then, the 
CER transition takes place once E exceeds a threshold 



TABLE I: Default values of parameters in the model. 

Here pm is a little larger than the experimental value (11 — 
14.5%) in Ref. [13] because the percolation threshold in 2D 
system is larger than that in 3D materials. The value of / 
is within the experimental range for manganites [3, Hlj]. The 
other parameters are relative values. 
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FIG. 2: (Color online) Resistances as a function of the applied 
field. Inset: &er as a function of /. 



Ec- In the subsequent stages of increasing/decreasing E', 
Ry remains low and nonlinear even for E well below E^. 
We find that the value of E^ increases with the IDI rate 
(not shown here), in agreement with the observations in 
Lao.225Pro.4Cao.375Mn03 [l3|, implying that the process 
is a relaxative evolution between metastable states. We 
use the ratio ^er = Ry^^'^ /Ry^^^ to characterize the 
CER, where Ry^^^ [Ry^^^) is the maximum (minimum) 
of Ry before (after) the CER transition. 5er appears to 
be proportional to /, as shown in the inset of Fig. 2. 

To understand the above transport behavior, four 
snapshots of PS pattern taken at different stages during 
the IDI process are shown in Fig. 3. Before E exceeds 
Ec, as shown in Fig. 3(a), the metallic sites accumulate 
into clusters in favor of lower interface energy. However, 
once E > Ec, the energy gained from dielectrophoresis 
overcomes the interface energy barrier, and the PS struc- 
ture is rapidly reconfigured into a stripe-like pattern, as 
shown in Fig. 3(b). As a result of the percolation, the 
system has higher overall conductance and higher effec- 
tive dielectric constant. Upon further increasing the field, 
the percolation path is strengthened, as shown in Fig. 
3(c), giving rise to a slight decrease of resistance. When 
the field is decreasing, the percolation path is hardly de- 
stroyed except when E is close to zero. When E is too 
small to compete against the surface tension and thermal 
energy, the percolation path will be corroded gradually 
and, given enough time, eventually break down, as shown 
in Fig. 3(d). Nevertheless, with the parameters as listed 
in Table I, the main shape of percolation path can re- 
main effective for a finite period of time, which accounts 
for the nonlinear transport behavior observed in Fig. 2 
after the voltage threshold is first exceeded. As a result, 
the dielectrophoresis in our simulation is irreversible. 

Besides the changing rate of the electric field, the in- 
terface energy coefficient A and temperature also play 
important roles in modulating the dielectrophoresis. In 
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FIG. 3; (Color Online) Four snapshots taken from an evolu- 
tion cycle of dielectrophoresis. 
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FIG. 4: (Color online) The field threshold as a function of 
the temperature. Here the surface tension coefficient A is 
chosen as the unit while the other parameters have the values 
in Table I. 



our simulation, the dependence of {Ec/A) upon (Td/A) 
is quite nontrivial, as shown in Fig. 4. {Ec/A) reaches 
the minimum at {Td/A) 1.1. It agrees with the ex- 
perimental evidences that there is a T-window for the 
dynamic PS (T ~ 25 — 85 K) 15| or a novel conso- 
lute critical point T ~ 30 K 13| in the phase-separated 
Lao.225Pro.4C!ao.375Mn03. The origin could be that, at 
large {Td/A), the thermal energy can destroy a stable 
percolation of the metallic phases; on the other hand, 
at small {Td/A), the relatively large surface tension can 
effectively prevent the metallic clusters from being re- 
shaped into stripes. Therefore, {Ec/A) is larger in both 
cases of large and small {Td/A). Although the parameter 
A can not be directly measured in experiments, the con- 
solute critical point does provide a way to compare our 
result of Ec with experimental results. Using the param- 
eters values taken from experiments 0, [l3 , 21 1 , in our 
model is estimated to be ~ 10^ V/m, within the same 
order of magnitude of the experimental data 13| . As we 
compare the physics of the present filamentous percola- 
tion picture and the CO-FM collapse picture, one notable 
thing is that Ec (~ 10** - 10^ V/m [lS[i3) required for 
the CO-FM collapse is indeed much larger than Ec in our 
model. The reason is that, the CER here is induced by 
merely repatterning of the CO and FM phases, not by 
converting them from one to the other, which costs more 
energy to compensate for the energy difference between 
CO and FM phases. 

Since the metallic filaments along the direction of elec- 
tric field break the space symmetry of PS, it is natural 
to expect transport anisotropy. In our simulation, re- 
sistances both along {Ry) and across {Rx) the electric 
field are affected, as shown in Fig. 5. Before the volt- 
age first reaches the threshold, the transport is isotropic. 



and resistances along both x and y are highly inde- 
pendent of the voltage drop. At the voltage threshold, 
the CER transition is accompanied by the sharp emer- 
gence of transport anisotropy. While Ry along the di- 
rection of the electric field is drastically reduced, a si- 
multaneous increase of Rx across the electric field occurs 
irreversibly near the threshold voltage. Furthermore, 
the system has the intriguing feature that it is metal- 
lic {dRy/dT > 0) along the electric field direction, but 
insulating {dRx/dT < 0) in the perpendicular direction. 

Discussion. Using the Bruggeman-Landauer (BL) 
equation based on self-consistent effective medium ap- 
proximation (SCEMA) [2^, the effective conductance of 
a well-distributed and isotropic binary mixtures system 
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FIG. 5: (Color online) Comparison of the perpendicular 
(cyan) and parallel (purple) resistances as a function of the 
applied field. 
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can be estimated as, 



1 



2(d- 1) 



{-P + ^P^ + ^{d-l)CiCM), (2) 



where d is the dimensionahty and C/ = \/Ri {Cm = 
^/Rm) is the conductance of insulting (metaUic) phases; 
P = CM{l-dpM) + Ci{dpM + l-d). Since Cj/Ca/ - 0, 
Eq. (2) can be simphfied to C w Ci/{1 — dpu) for 
all cases (/ ^ 1) in our model. Therefore, the ef- 
fective resistance is about (1 — dpM)Ri — 0-6Ri for 
all values of / considered here and the conductance is 
isotropic. In our simulation, the conductance is isotropic 
and its value is indeed independent of / before the 
CER transition. The value of Ry w 0.55i?/ is close 
to, but slightly smaller than, the BL estimation. The 
reason is that the accumulation of metallic phases into 
clusters will enhance the conductance. As the perco- 
lation path is formed, the Ry^^^ can be estimated as 
RiKpm! + 1 - Pm) ~ Ri/ipAif) since / > 1. It ac- 
counts for the linear dependence of 5er on /, as shown 
in the inset of Fig. 2. (The relationship Seb. oc Ri/Rm is 
independent of the choice of Bm /^i in our model.) On the 
other hand, R^'^^ after the CER transition is approxi- 
mately — Pm) + RmPm ~ Ri{^ — Pm)- Compared 
with its original value, there is a relative small change 
(increase) of R,: = (1 -m)/(l - 2m) 

by the two-dimensional BL equation. 

Compared with the dielectrophoresis in liquids, there 
are two remarkable features in phase-separated mangan- 
ites, adding to its appeal to the research community. 
Firstly, the prototype PS in manganites would gener- 
ate inhomogeneous electric field by itself, while for liq- 
uids great effort was made to artificially produce non- 
uniform electric field [1, Yet, as in the cases of the 
liquids, the PS in manganites can also be controlled by 
applying artificially non-uniform electric field (e.g., pat- 
terning of gates), offering great potential for making de- 
vices. Secondly, the diversity of electric properties be- 
tween different phases in manganites, e.g. resistance and 
dielectric constant, is far greater than those in conso- 
lute liquids. In real manganites, the colossal electroresis- 
tance/magnetoresistance implies a colossal difference be- 
tween the resistivities of metallic and insulating phases. 
Therefore, the main difhculties of dielectrophoresis in liq- 
uid, i.e. a large voltage required and weak effects, could 
be overcome in manganites. 

In conclusion, we introduced the dielectrophoresis dy- 
namics, a phenomenon usually existing in mixed phase 
liquids, to investigate the phase separation in solid man- 
ganites. We have shown that the colossal electroresis- 
tance could be achieved by the filamentous percolation, 
through a dielectrophoresis dynamic process without in- 
creasing the metallic composition. Our results give a 
probable explanation of the colossal electroresistance in 
manganites, including its time-dependence, irreversibil- 



ity and the consolute critical point. Besides, it was also 
predicted that the dielectrophoresis process renders the 
system anisotropic. All these properties can not be ex- 
plained by the physics of single phases but arise from 
phase competition. The present work also highlights the 
analogy [sj between strongly correlated electronic sys- 
tems and complex systems in the classical world, e.g., 
complex fluids. 
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